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A theory is presented for calculating the effect of the electromagnetic field on the centre of mass 
of a macroscopic dielectric body that is valid in both quantum and classical regimes. We apply the 
theory to find the classical equation of motion for the centre of mass of a macroscopic object in 
a classical field, and the spreading of an initially localized wave-packet representing the centre of 
mass of a small object, in a quantum field. The classical force is found to be consistent with the 
identification of the Abraham momentum with the mechanical momentum of light, and the motion 
of the wave-packet is found to be subject to an acceleration due to the Casimir force, and a time 
dependent fluctuating motion due the creation of pairs of excitations within the object. The theory 
is valid for any dielectric that has susceptibilities satisfying the Kramers-Kronig relations, and is not 
subject to arguments regarding the form of the electromagnetic energy-momentum tensor within a 
medium. 
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In recent years there has been significant progress in 
cooling macroscopic degrees of freedom down to temper- 
atures where deviations from classical behaviour can be 
observed pM]. This is to the extent that when suffi- 
ciently isolated from other degrees of freedom, the macro- 
scopic mechanical motion can be reduced to the quantum 
mechanical ground state [5HZ]- It thus becomes feasible 
to prepare superpositions of macroscopic states of mat- 
ter [SI [9] , where quantum mechanics is known to be at 
its most perplexing, and where the wave-function is par- 
ticularly hard to take as a serious description of real- 
ity maim. 

Motivated by this remarkable work we look to con- 
tribute through asking how we should properly describe 
the motion of a macroscopic object interacting with a 
quantum field. To be more specific, suppose that we 
have an object composed of ~ 10 12 atoms, interacting 
with one another and the outside world through the elec- 
tromagnetic field, and where the object may have been 
prepared in such a way that its centre of mass, R can 
be considered as a quantum mechanical variable with an 
associated operator, R (see e.g. P21 E]). We wish to 
understand the dynamics of the centre of mass in this 
situation. 

It is evident that we cannot make any progress if we 
work in terms of a microscopic description. Yet, the com- 
paratively macroscopic scale of the object has the ad- 
vantage that we can assume that the collective internal 
atomic motion is close to some equilibrium configuration, 
and that any linear perturbation of the system away from 
this can be characterized via a susceptibility. We there- 
fore look for a Hamiltonian that can be derived from an 
action principle, correctly describes the dynamics of the 
centre of mass and the electromagnetic field, and is for- 
mulated in terms of susceptibilities rather than atoms. 
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This is certainly not the first attempt to tackle this sort 
of problem. Some time ago C. K. Law worked out the 
theory of the quantized electromagnetic field in a cavity 
when one of the walls moves according to the rules of 
quantum mechanics [14] . More recently there have been 
several papers investigating the quantum motion of the 
centre of mass of macroscopic dielectric objects within 
the electromagnetic field [T5hT7] . often treating the sys- 
tem in the context of the formalism of open quantum 
systems, and considering the interaction of the object 
with a single cavity mode. 

We emphasize that we do not yet aim to find a good 
model for the cavity configuration, but rather look for 
a theory that correctly treats the motion of macroscopic 
bodies in an electromagnetic field. One major new in- 
gredient we add to previous work is that we develop a 
description that fully includes the effects of dispersion 
and dissipation of the electromagnetic field within the 
body. As a byproduct we shall find that we can also 
provide unambiguous expressions for classical forces on 
dielectric bodies, with results that may well be relevant 
to the continuing discussion of optical momentum in me- 
dia 18 -25j. In this context our work is perhaps most 
closely related to |19) . 

One may alternatively wish to see the following work 
as a kind of generalization of Casimir-Lifshitz the- 
ory [211 [27] to the situation where the centre of mass 
of the dielectric is treated as a dynamical variable. Our 
findings may clarify the difference between the results 
of [28] and [2Z[[29]. 



I. ACTION AND CLASSICAL MOTION 

We begin with the action necessary to describe the 
classical motion of a dielectric in an electromagnetic field. 
The Lagrangian associated with the free non-relativistic 
motion of the centre of mass of a body must take the 
form, L cu = ^M(dR/dt) 2 , where M is its mass, and 
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dR/dt is the centre of mass velocity. Similarly the La- 
grangian associated with the free electromagnetic field is 
also well known to be, Jzf F = (e /2)(E 2 -c?B 2 ) [3D]. The 
challenge is in how to characterize their interaction. 

The key to understanding the interaction between field 
and dielectric can be found in the work of Huttncr 
and Barnett |31j . and the later more general findings of 
Philbin [3H [33J . They find that the dynamics of a sys- 
tem that is characterised by a linear susceptibility can be 
precisely mimicked using a field of simple harmonic oscil- 
lators. To understand why this must be so, consider the 
polarization field of a stationary dielectric that responds 
linearly to the electric field, 



P(x,t) 
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Due to the fact that P is a function of the electric field in 
the past and not the future, x EE (f2) must be an analytic 
function of frequency in the upper half of the complex 
plane [3U [35]. This has the consequence that the real 
and imaginary parts of x EE & r e related by the Kramers- 
Kronig relations 
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where 'P' indicates the principal part of the integral. In- 
serting ([2| into (JlJ and using the convolution theorem, 
we find that the polarization field can be written as 
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P(x,t) = dw 
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Within pi we have the quantity, G u (t - t') = 8(i - 
t')sm[ijj(t — t')]/oj, which is recognised as the retarded 
Green function for a simple harmonic oscillator of unit 
mass and natural frequency w [HI [36] . Expression ([3| 
could also be obtained if the polarization field were iden- 
tified with an integral over a continuum of simple har- 
monic oscillators, all coupled linearly to the electric field, 
each oscillator having a mass (2wIm[x E E(w)]/7r) _1 . 

With this in mind, the Lagrangian density describing 
a stationary dielectric in an electromagnetic field is 



^ = ^ F + P-E + M-B 
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where A takes two indices, E and B, these labelling the 
two sets of oscillators [32] . The polarization and magne- 
tization are given by 
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where a(u>) — 
v /2wIm[xBB(w)]/7T. 



A/2wIm[x EE (w)]/7r and /3(lj) = 
Note the equivalence between the 



imaginary parts of the susceptibilities appearing within 
the mass of the oscillator field (as in the previous para- 
graph), and appearing within the coupling, as in 
From this action a Hamiltonian can be derived, and the 
system can be canonically quantized. Recent work has 
found Q to provide a rigorous basis for the theory of the 
Casimir effect [32] 

When the dielectric is set in motion, the polarization 
and magnetization are no longer given by ([5]), but to 
first order in R transform as, P 1 = P + R X M/c 2 , 
and M' = M Rx P. Furthermore, the susceptibility 
becomes non-local, and a function of the Doppler shifted 
frequency, f2' = fi — R • k [35] . Taking these factors into 
account, the action for a uniformly moving dielectric was 
derived in [37] and from this a quantum theory moving 
media was constructed 38 . The action was found to be 
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where we write the electromagnetic field in terms of the 
four-vector potential, E = — Vtp — A, B = ~V X A. The 
Lagrangian density associated with the two sets of simple 
harmonic oscillators is modified from the stationary case, 
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and the interaction Lagrangian density is also changed, 
due to the aforementioned transformation of the polar- 
ization and the magnetization, 

POO 

.s?int = / duj \- E • + B ■ «baM]-xL a) (8) 

The a X \' matrices being given by, a EE = l3a(w); a BB = 
i 3j 5(w); « EB = (l/c 2 )-R X ± 3 /3(u); and a BE = R X 
l3a(w). In the above expressions we have suppressed the 
spatial dependence of the susceptibilities, which depend 
on the difference between the integration variable, x and 
the position of the centre of mass, R, i.e. x EE = X EE ( a3 — 
R,uj), being zero everywhere outside the dielectric. 

A comparison of expressions (|6j^8j) and those of [33 [3S] 
will show that we are working under the approximation 
that (R/c) 2 ~ 0. This is because we will now treat the 
motion of the centre of mass as a dynamical variable 
rather than an external parameter. Beyond first order in 
R/c the dynamics of the centre of mass of an object are 
extremely subtle. Properly speaking the centre of mass 
becomes a redundant variable, and one must instead deal 
with the centre of energy [30l I39H42] . 

As a final comment we note that we only consider the 
translational motion of the body, and not its rotation 
about the centre. It is possible to include this through 
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introducing the three Euler angles into spatial depen- 
dence of the quantities, ol\\i , but we do not study these 
effects here. 



tensorial quantity playing the role of the mass, fj, con- 
tains contributions from both the intertial mass of the 
dielectric and the spatial distribution of the oscillator 
amplitudes, 



A. Electromagnetic force on a dielectric in free 
space 

The action given in ^ allows us to calculate an un- 
ambiguous expression for the force on the centre of mass 
of a dielectric body. To do this we need only add the 
kinetic energy term, L CM to the Lagrangian, 



.s" = / L'dt = s + \ I MR 2 dt 
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Then the classical equations of motion of the centre of 
mass can be obtained from an application of the Eulcr- 
Lagrange equations to the variable, R(t) 



d (dL'\ _ dV 
dt \dRj ~ dR 
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We find that the associated canonical momentum is given 
by 



li = Ml 3 + 
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du;(V®Xy).(Xy®5) 
(12) 

The second term in ( |11| — involving a coupling of the os- 
cillators and the fields — was also found in However, 



the second line of (11) and second term within fi ap- 



pear to be new. The origin of these terms can be traced 
back to the continuum description of the medium, and 
ultimately the Doppler shift within the Kramer-Kronig 
relations. For example such a contribution is not present 
within the single oscillator model of the internal degrees 
of freedom given in |43| . 



To calculate the right hand side of ( 10 ) we note that 



the Lagrangian depends on R only through the coupling 
matrices, a\\> within L INT = J Jzf INT d 3 £c. The derivative 
of this quantity with respect to the zth component of the 
centre of mass vector is found to be 



dV . 

p — — — ii- R 
dR 
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which includes a contribution from the both the harmonic 
oscillators alone, and the coupling of the electromagnetic 
field to the oscillators. The polarization and magnetiza- 
tion are the rest frame expressions given in ([5]), and the 
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(13) 

Due to the fact that the coupling only depends on the dif- 
ference between the integration variable, x and the centre 
of mass position we can equivalently write, dct\\> /dRi = 
— dct\\i/dXi. Inserting (111 and (13) into (10) then gives 
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an expression for the electromagnetic force on a moving 
dielectric body, 
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To obtain |14[ ), which does not make explicit reference 
to the oscillator fields, we integrated out a term equal 
to the gradient of a scalar containing only oscillator am- 
plitudes and not the electromagnetic field. In the region 
outside the body where the a.\\' = 0, the oscillator field 
is assumed to be at rest, therefore the divergence of any 
quantity containing them will integrate to zero. We also 
applied the equations of motion for the two sets of oscil- 



lators, 

d_ 
dt 



RV 



=ctl x -E + cxl x -B (15) 



To further simplify ( 14 ) we need to impose the equa- 



tions of motion for the electromagnetic field. These fol- 
low from a variation of the action ^ with respect to the 
four-potential, = ((p/c,A), and are the familiar ex- 
pressions; V-D = 0;V'B = 0;Vx£ = -dB/dt; and 
V X H = D. The auxiliary fields being given by D = 
e E+P+(l/c 2 )RxM andH = ^B-M+RxP [55]. 
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Applying these equations of motion to ( 14 1 gives 
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where V is the volume of the dielectric, dV is the two 
dimensional region of free space on the boundary of V , 
and dS is an infinitesimal element of surface area directed 
outwards from the centre of the body. 

In the case wh ere the body is instantaneously at rest, 
R = 0, then ( 16 ) can be reduced to the compact form, 



MR- 
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where the momentum density is given by the Abraham 
expression [20] . 



V A = ~E X H 



and the stress tensor is of the Maxwell form 



(18) 
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B 



- 2 B ■ B) 



(19) 

In deriving (17) we have neglected a term proportional 



to the acceleration of the body times its spatial extent 
divided by c 2 . This is required within the present con- 
text, for the acceleration must be small enough such that 
we do not have to worry about relativistic effects e.g. (44j 
Chapter 6]. 

The expression for the acceleration of the centre of 
mass ( [17] ) appears to agree with existing findings. Con- 
sider the case of a short pulse entering an inhomogeneous 
dielectric. During the time it propagates through the 
body, the field outside is approximately zero and expres- 



sion (17) suggests that the centre of mass should ex- 



perience a force equal to minus the rate of change of 
the Abraham momentum within the medium. This is in 
agreement with the result that the mechanical momen- 
tum of the pulse within the medium is associated with 
the Abraham expression [23, [24]. Furthermore, if we con- 
sider an optical field interacting with the medium, and 
time average (17) over an interval much longer than an 



optical cycle, then the force simply equals the integral of 
the averaged Maxwell stress tensor over the surface of the 
dielectric, which is the familiar expression used in studies 
of radiation pressure in free space [IS], and is equivalent 
to the average net macroscopic Lorentz force on the body 
(see [46]). 



Although (17) is in itself not a surprising result, it is 



the procedure of constructing an energy momentum ten- 
sor as is done in e.g. [151 H3- Instead we started from 
an action that was constructed to give the correct con- 
stitutive relations for a moving medium and then varied 
this action with respect to the centre of mass velocity. 
We might understand this procedure as being similar to 
the likes of [47] 08] , where the Doppler shift in the con- 
stitutive relations of a medium is shown to be intimately 
connected with the force imparted by radiation. This 
approach appears to be free from the ambiguity associ- 
ated with the usual division of the full energy momentum 
tensor into 'field' and 'matter' contributions [19] . 

If the reader performs the necessary manipulations to 



pass from (14) to (17) they may notice the divergence of 



interesting that it can be derived without carrying out 



the Minkowskii stress tensor, er M = D(g)E + B§§H — 
(1/2)1 3 (£) • E + B • H) appears within the intermedi- 
ate expressions (before divergences of terms containing 
the polarization and magnetization are dropped). It is 
tempting to conclude that this tensor is somehow rele- 
vant to the force density within a medium. However such 
a conclusion is not at all justified in the case of an isolated 
body in free space, given that once integrated over the 
volume of the medium there is no meaningful distinction 
between this intermediate tensor and ( |l9| ). To investi- 
gate this point further, and to illustrate the application 
of our formalism to more than one dielectric body we 
now examine the electromagnetic force at the interface 
between two media. 



B. Electromagnetic force at an interface between 
two dielectric media 



Consider the specific configuration shown in figure 
[l] where we have two planar dielectric media instanta- 
neously at rest and in contact over the x — plane. The 
first medium occupies the half space x < 0, and is taken 
to be of a mass so great that we can consider its posi- 
tion to be fixed. The second medium occupies the space 
< x < a and is allowed to move in response to its inter- 
action with the electromagnetic field, which will include 
some contribution from the field at x = 0. 

To describe this situation in terms of an action we must 
introduce another pair of oscillator amplitudes into d9]). 
We denote these by , and they describe the collective 
degrees of freedom of the medium occupying x < 0. The 
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FIG. 1: Two planar dielectric media are in contact over the 
plane x = and the centre of mass of the first medium, R(a > 
x > 0) is free to move in response to the electromagnetic field. 
The internal degrees of freedom of the medium occupying 
a > x > are represented by the two sets of oscillators, 
, and those of the medium occupying x < by y£ A ' . 



action given in ^ is then equal to, 

S" = S> + J d+x J™ do,{ £ I (y[ X)2 - u?Y 
+ C(u)E ■ Y <f) + fl ■ yW | - J V(R)dt (20) 



r(\)2 



Both media are assumed again isotropic and to be ade- 
quately represented by e and /i in the rest frame. The 
coupling constants, C( w ) an d perform the same role 
as the a(uj) and /3(w) in |5l), but are non-zero in the re- 
gion x < rather than < a; < a. The potential function 
V(-R) represents the short range microscopic fields that 
prevent the two media from overlapping. One obvious 
choice for V{R) would be an infinite potential barrier in 
the region x ■ R < a/2. 

As R and R do not appear within the new terms in 



(201, the equation of motion (14) is unchanged except 
for the addition of a term equal to —'VrV(R) on the 
right hand side. However, we should be careful not to 
make the identifications P and M within (fl4| as given 



by ([5]). The total polarization and magnetization that 
appear with Maxwell's equations are now given by, 



M = 



a( W )xW + CM^ E) 
P(u)XW+Z(u)YW 



duj = P 1 + P 2 

duj = Mi + M 2 

(21) 



manner as in the previous section we are left with 
MR+~ I Ex{B/iMi-M x )^x + V R V{R) 



V • a + (V • P 2 )E - (P 2 + V X M 2 J X B 

(22) 

To proceed we note that the following identity can be 
derived from Maxwell's equations, 



d 



V • a - e„— (E X B) = -(V • P)E 



dt 



P + V X M x B (23) 



After applying (23) to (22) we find that the force on the 



centre of mass of the dielectric again obeys the equation 
of motion given by (17), plus the short-ranged force due 
to the potential, V(R). In this case too the theory pre- 
dicts the centre of mass of the body to accelerate accord- 
ing to the net macroscopic Lorentz force, with the Abra- 
ham momentum relevant for the mechanical momentum 
exchanged with the electromagnetic field. To prevent 
misunderstanding we should emphasise that this method 
of including dispersion and dissipation would not be valid 
for a body surrounded by a fluid - nor for the deformation 
of an elastic body - as the relevant action should include 
the local degrees of freedom of the body (e.g. the strain, 
or the fluid velocity), which we do not. Therefore this 
finding is not immediately at odds with [23 IH] . Nev- 
ertheless, in the remaining work we confine ourselves to 
the uncontroversial regime where the media are not in 
contact, and are surrounded by vacuum. 



II. QUANTUM THEORY OF MOTION 

Having established that the action ^ reproduces the 
known results of the classical theory of radiation pres- 
sure, and mechanical momentum of light within media, 
we construct the corresponding quantum theory. For the 
sake of simplicity we construct the Hamiltonian for a sin- 
gle dielectric body isolated in vacuum. The case of two 
or more interacting bodies is not fundamentally different, 
and one can proceed as in the previous section. 

As usual, first find the canonical variables. The canon- 
ical momentum for the electromagnetic field is given by, 



dC 
dA 



-e E - P 



1 
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and for the oscillator fields, 



n 



(A) _ 
X u 



dX 



dC 
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(24) 



(25) 



we should therefore amend ( 14 ) with P — > Pi and M — > 
Mi. After applying Maxwell's equations in the same 



The canonical momentum associated with the centre of 
mass coordinate is given in (11). Having obtained these 
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variables, the Hamiltonian is then constricted in the 
usual way, 



H = J d 3 x (n A + dojll^l ■ X 



+ P-R-L' (26) 

However, the relationship between the canonical vari- 
ables and the mechanical ones is a little involved. To 
cast the Hamiltonian entirely in terms of canonical vari- 
ables we use the relation, 



la 



Mo 
Mc 2 



d 3 x (M 2 1 3 - M ® M) 



■A = ^(P + r) 

(27) 



where, 



r = J d 3 x[n Q M x (n A + P) + p x (v x A) 

-^|°°^(v®xW).nW] (28) 



To simplify (27), we might assume that the energy equal 



to the integral of the square of the magnetization over 
the medium times fio is vanishingly small in comparison 
to the rest energy, Mc 2 . From the perspective of classical 
physics this is permissible, but quantum mechanically the 
magnetization energy will contain a ground state contri- 
bution, which may in total not be small. To account 
for this possibility, we assume an isotropic body and 
absorb this approximately constant energy — whatever it 
is — into the rest mass, M of the body. Then (27) be- 
comes, R = (p + r)/M, and (26) is given in terms of 
canonical variables by, 



H = 



(p + T) 



2M 
M-B + 



d 3 x 



Zeo 



It, 



dcu n 



-uj a X 



(29) 



where we have applied the constraint that V • 11,4 = 0, 
and again dropped terms of order (R/c) 2 . Comparing 
(29) to the Hamiltonian of [32] it is evident that the in- 



tegral to the right of the kinetic energy term is the ex- 
pression for the Hamiltonian of the electromagnetic field 
interacting with a dielectric at rest. The centre of mass 
motion is coupled to the field and the dielectric by the 
term T, given in (28). This term is analogous to the 



vector potential in the Hamiltonian of a charged point 
particle interacting with the electromagnetic field. How- 
ever this 'vector potential' is given by a quite complicated 
expression involving an integral of both field and material 
degrees of freedom over the whole body. 

Comparing p9| with the expressions of Q31 [T5] we see 
that Law and Cheung's result — that a dielectric body 
couples 'minimally' to the electromagnetic field — also 
holds when dispersion and dissipation are accounted for. 
Given this form of the coupling, we expect Aharonov- 
Bohm type effects to be present within the behaviour 



of the centre of mass, with the semiclassical phase shift 
given by, 



A0 = - 



TdR 



(30) 



Presumably this phase is very difficult to directly observe 
within an interference pattern — e.g. in the same manner 
as for neutrons [50] — and some thought must go into the 
indirect implications. We note that the analogue of this 
effect was also found in [16] for the case where dissipation 
and dispersion are neglected. The form of the effective 
vector potential given in ( 28 1 contains terms that may be 
identified as the equivalent of the Aharonov-Casher [51] 
and He-McKellar-Wilkens [52l [53] contributions for an 
extended object. However, we find an additional term 
is also present due to our continuum description of dis- 
persion and dissipation. This does not arise in the case 
where the object is treated as a discrete collection of os- 
cillators [33], but seems to be necessary in macroscopic 
electromagnetism. 

The Hamiltonian given in ( 29 1 is bounded from be- 
low. In |38j it was found that if we attempt to construct 
a quantum field theory for electromagnetism interacting 
with a uniformly moving dielectric then the Hamiltonian 
is extremely odd and, without including the effects of spa- 
tial dispersion, unbounded from below. We again inter- 
pret these odd features as arising from a failure to include 
the dynamics of the centre of mass motion, which are in- 
separable from those of the field. The Hamiltonian of [38] 
should therefore be understood as an approximate de- 
scription of the physical situation, with ( 29 1 being more 
fundamental. 

The quantum mechanical Hamiltonian is obtained 
from p9l in the usual way. The canonical variables be- 



come operators satisfying the canonical commutation re- 
lations. The centre of mass operators obey, 



R,P 



(31) 



and the electromagnetic field operators satisfy, 

A(x,t),fL A (x',tj\ = ih8 A _(x-x') (32) 

where S±(x — x') is the transverse delta function [54] . 
Finally, the commutation relations for the oscillator field 
operators are, 



X { *\x,i),TL { x\{x',t) 



= iht 3 5 xw 5{uj-J)8 {3 \x-x') (33) 
The Hamiltonian operator then takes a very similar form 



to (291 
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where, 
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(35) 



the only distinction being that a symmetric ordering of 
the oscillator field operators and their associated canon- 



ical momenta must be chosen in the T operator, 



d 3 xi fi M X 



(p-A 



P x (V x A) 
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(X 



(36) 



The part of the (34) given by (35) is identical to that 'polaritons'. 



of a stationary dielectric at a position given by R, and 
may be 'diagonalized' into 'normal modes' ('polaritons' 
in the sense of [3T]) of the combined system of the field 
and the internal degrees of freedom of the dielectric [35] . 
The total 'polariton' energy can then be written in 
terms of a number operator. The structure of (34) is 



made especially interesting by this fact that the number 
operator depends on the centre of mass operator, R: for 
each possible position of the dielectric body we have a 
different basis of 'polariton' number states. 

From this discussion, it is not clear what the ground 
state of (34) is. We might expect that the ground 



state would correspond to a completely de-localized 
centre of mass coordinate and zero 'polaritons', i.e. 
IV-'o) = ^Vo|0)hj where Nq is a normalization factor, and 
|0)h is the state of zero 'polaritons' when counted using 
the basis associated with the centre of mass being at J?. 
However, this expression for does not appear to be 
an eigenstate of (34), due to the fact that p\0)r ^ 0. We 



leave the problem of finding these eigenstates for future 



work, and to conclude we apply ( 34 ) to the spreading of 
a localised centre of mass when we initially have zero 



Quantum motion of an initially localised centre 
of mass 



In order to glean some understanding of the quantum 
centre of mass motion, we consider a state where at a time 
t = the centre of mass has been prepared so that it is 
sharply localised around some mean position, R (see fig- 
ure^. If the mass of the object is large enough then the 
characteristic spreading time of a wavepacket of width 
a -1 / 2 , T s — M /ha can be long compared to the dynam- 
ics of the field. For example, were we to know the initial 
centre of mass position to within 10 _9 m (a ~ 10 18 m 2 ), 
then a system of ~ 10 12 atoms (M ~ 10 -14 kg) would ex- 
hibit T s ~ 10 2 s. Therefore, in the regime of macroscopic 
electromagnetism, a good approximation is to expand 
(34) around Rq [T3], a significant time being required 
before terms of order (R — Rq) 2 and higher become rel- 
evant to the time evolution of \tp) [56] . Such a first order 



expansion of ( 34 ) yields 



H 



P 



2M 



(R -R )- J d 3 x ^(V fl (2) P) • (TLa + Po) - (Vr 8 M)o • (v X A) 



H (37) 



where the zero subscript indicates that the quantity is 
evaluated at R = Rq, Hq = H p (Rq), and the ellipsis 
denote further terms in the expansion of T around Rq. 
As mentioned in the previous section, we can diagonalise 
Hq into 'polariton' modes with the transformation of |32j 



(see appendix [A]) , which gives 



+ h.c. (38) 
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We shall often assume that Hq takes the above diago- 
nalised form, and that the remaining 'polariton' (non- 
mechanical) degrees of freedom in (37) are expanded in 
terms of the C u operators. 

Before proceeding with the calculation, we choose to 
transform (37) into a different representation. Due to the 



independence of ( 38 1 from R we can quite simply cast 



(37) into an interaction representation via the unitary 



transformation, 

IJ\ (t) = cxp 

where we identify, 

\1>) = &1\<P) 



it 



H = Ut (Hz + Ho) U\ 



(39) 



We then find the usual result, Hi\ip) = ihd\ip)/dt, where 
Hi = U\{H — H )U~i, and the effect of the unitary trans- 



formation on (37) is to make the electromagnetic field, 



polarization, and magnetization operators all time de- 
pendent, via the substitution C u —> C u exp (— iujt). In 
order to simplify the situation further, we make a second 
unitary transformation, 



U 2 {R,t) =exp 



-j_(R-Ro) 



(40) 



and Hi = U 2 H' I U\ 
representation of the 



where we write, \<p) = |x) , 
ih(dU 2 /dt)Ul This is the 
Schrodinger equation which we apply to the motion of 
the wave-packet: H'j\x) — ihd\x)/dt, the approximate 
Hamiltonian being given by, 



H'i(t) = ^ + (R-Ro)-l / d 3 x 



1 



(V H ® P)o • (Ha + Po) - (V fi ® M) ■ V X A 



1 

2M 




P ■ (T ® Vjj)o • (R - Ro) + (R— Ro) ■ (V fl ® f ) • pi (41) 



where all terms have been truncated to first order in 
R — i?o, and the Hcrmitian part of the resulting ex- 
pression has been taken. When evaluating the average 
acceleration of the wave-packet, the second line of (41) 



will in general give something distinct from the corre- 
sponding classical result. This is because the ordering of 
the position and momentum operators matters. However 
we do not investigate such effects here. 

Our attention now turns more precisely the dynamics 
of a wave-packet initially prepared in the state illustrated 
in figure [2] 



|X(* = 0)>=(£) 



3/4 



r ( R— Ro ) z 
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(42) 



where |C) is an eigenstate of the C w (x, Rq) operators. 
It is reasonably straightforward to show that for such a 



state, the average position is subject to a mean instan- 
taneous acceleration that can be interpreted as the ana- 
logue of the classical result (17). Firstly, via the relation 



(R) = —(i/h)([R, H'j]) wc find the average velocity of 
the state (42) to be zero, 



R) = 4 <xl P + (R - Ro) ■ (Vr (2) f ) 



AI 



lx) = o 

(43) 

Introducing the quantity, V = A/ -1 [p+(R— Rq) ■ (Vjt® 
r)o], the average mechanical force can be similarly com- 
puted from the relation, M(R) = -(iM/H)([V, H'j]), 
which yields, 




R ' 



tVq x M } • E 



V R ®P) •E+(V R ®M) -B 



M - Vn X P • B 



dt 



dT 

o dt 



\x) 



v-(v H «r)o |x) 



(44) 



where we have defined the quantities, E = -e Q + an d applied the result (|43|. The velocity dependent con- 

P +cr 2 VxM Q ], E = -e-^n^+Po], and-B = VxA, 
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FIG. 2: A small object, characterized by e and /j, satisfying 
the Kramers-Kronig relations is prepared so that the proba- 
bility densit y fo r the centre of mass is spread over a width, 



a ' (c.f. (42 1). The coupled system of field and dielectric 



is initially in the state \C), and we investigate the effect of the 
coupling between this system and the centre of mass on the 
time evolution of the wave-packet. 



tribution to the electric field operator drops out in the 
second line of (44), due to the form of (|42|). Note that the 



expression on the first line of (44) is the operator equiv- 
alent of (10). Besides assuming that the field does not 



vary significantly over the extent of the wave-packet, the 
derivation of ( 44 ) only relies on the fact that ([42 ) is both 
separable and rotationally symmetric around the point 
Ro . Applying the equations of motion of the operato rs— 
dTo/dt = -(i/H)[r ,H ]— to the time derivative in Q , 



and performing several integrations by parts, one obtains 
the result, 



M(R) = ( X \ 



where "P / 



<r • dS — 

dv c 



V A d 3 x)\ X ) (45) 



E X H, H = B/hq - M , and <r 



-2, 



e Q [E®E+c 2 B®B-{l/2)t 3 (E +c 2 B )]. It is therefore 
clear that the average position of an initially stationary 
and well localised centre of mass is subject to the expec- 
tation value of the operator equivalent of (17) [57 . This 
is comforting from the point of verifying that this theory 
gives sensible answers, but is wholly anticipated given 



that expectation values are known to exhibit behaviour 
consistent with classical equations of motion. 

We now move from a discussion of average properties 
computed from the state, to the state itself. To solve the 
Schrodinger equation we proceed pertubatively, the state 
being expanded as \x(t)) — lx (*)) + Ix (*)}) obeying 
the zcroth and first order relations, 



P 

2M 

£. 

2M 



x (0) W> 



\x (1) (t))-m 



ih 



d\x W (t)) 

at 

: ,d\X (1 Ht)) 



dt 



= -AH(t)\x W (t)) (46) 



where H'j(t) = p 2 /2M + AH(t). The idea is that, as 
alluded to at the beginning of this section, if the spread- 
ing time of the wave-packet is long then over some rel- 
atively extended time interval the object will remain 
well localised around Rq. With some approximation we 
can then consider the first power of R — Rq to domi- 
nate the time evolution of \x)- Applyi ng the condition, 

|X(* = 
tion, 



0)) = \x (0) {t = 0)), the first of (46) has the solu- 



\x {0) (t)) = 



1 



•h ihctt/M J 



3/2 



g 2(l + ihct/M) \Q 



(47) 



After inserting ( 47 ) into the second line of ( 46 ) we can use 
the free particle Green function, K(R, R'\T){T = t-t') 
to find an expression for |x^(^))- The free particle Green 
function is a solution of, 



E__. h 9_ 
2M dt 



K(R, R'\T) = -iM (3) (R - R')S(t - t') 



and has the explicit form, 

K(R,R';T) = (2^) 3/2 ^ (R ^' )2 ©(r) (48) 



which is the retarded Green function. In terms of (48), 
the solution to the second line of (46) is, 



lx (1) W) = 



J dt 1 J d 3 R'K(R,R';T) 

xAH(R',1/)\x^(R'^)) (49) 



The spatial integral in (49) is Gaussian, which can be 



evaluated using the usual results. Inserting (]41h , (47) 



and (48) into (49), and performing this integral gives, 



x (1) W) = 



dt' 



av, 

dt' 



dS-. 



av 



X (0) W» (50) 
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where we have applied the field operator equations of tations as the wave-packet spreads (i.e. there are pairs of 



motion in the same way as in the derivation of ( 45 1 , and 
kept to first order in R — R . We note the appearance of 
(Vfl-r)o within the curly brackets in ( 50 ). Such a term is 
certainly peculiar to quantum mechanics, where it arises 
because position and momentum variables do not com- 
mute. This term stems from the fact that the canonical 
momentum of the centre of mass has both a contribution 
from the velocity of the medium and the field degrees of 
freedom. The second term in curly brackets is the op- 
erator equivalent of the classical radiation pressure force 
investigated in section \L~K\ Being multiplied by R — Rq, 
this serves to shift the centre of the wave-packet in accor- 
dance with (45). However, besides shifting the average 
position of the wave-packet and leaving the field in the 
|C) state, this term also creates pairs of 'polariton' exci- 



'polariton' creation operators within (50)), and this will 



cause a fluctuation of the form of the packet, as well as 
a change in the state of the field. As a final calculation 
we look to make this discussion more concrete and dis- 
cern the regime of parameters where the above effect is 
evident. 



The field operators in ( 50 ) are expanded in terms of 
the Coj and Cj, operators given in appendix [a] We also 
specialise to the case where \Q = |0), which is the state 
for which C w (a;)|0) = 0, for any choice of annihilation 
operator. Applying ( |A16 - A17) we find that the operator 
equivalent of the stress tensor acting on the vacuum state 
gives a sum of two identifiable terms, 



poo pOO f -| ^ 

a(x)\0)=eoj dSl J dO'e i ( n+n '^' li)(x, SI, SI') - -l 3 Tr [v(x, SI, fi')]j |0) 

+ — lim / dnim<T(x,x'n) - il 3 Tr[T(x,x',fi)]l |0) (51) 

7T x^x' J Q I 2 J 



where in ( |51[ ) we have defined 

t(x, x',si) = n 2 G(x. x'\ si) + c 2 v x g(x, x'-, si) x V' 

(52) 

and 

v(x, SI, SI') = e\x, SI) ® s\x, SI') 

+ c 2 b\x,SI) ®b\x,S1') (53) 



We can interpret (51) in a reasonably straightforward 



manner. The final time-independent term is identical to 
the vacuum stress tensor used in Casimir-force calcula- 
tions (see e.g. (551 HZ]). This term is multiplied by R—Rq 



in (50), and consequently shifts the peak of the wave- 
packet in accordance with the Casimir force. Clearly this 
force would be expected to vanish in the case of a lone 
body in free space. Meanwhile, the first line is associated 
with the production of pairs of excitations propagating 
from the surface to within the volume of the dielectric 
body (c.f. ( [A15] )) 
body in vacuum. 



and is non-zero even for an isolated 



which is a well known feature of the vacuum stress ten- 
sor [26]. Yet (50) contains the integral of <r|0) over a 



closed surface, and the divergent contribution to (51 ) of- 
ten vanishes when integrated over a closed surface. The 
standard procedure for dealing with the infinite contri- 
bution is to subtract a Green function for a homogeneous 
medium from the Green functions in (52). This proce- 
dure is justified in appendix [5] where it is also shown 
that the vacuum force on a lone object vanishes in free 
space. Hence the second line of (51 1 does not contribute 
to ((50). 



In addition to a, (50) also contains the time derivative 
of Va\0), which is, 



or 



oo poo 



JO 



*\0)=i dSl dSl'e l< - n+n,)t '(Sl + Sl') 



As it stands, the term on the second line of (51) is 
not finite. This is because the integral over SI diverges, tion over time, we find, 



£ \x,Sl) x n^{x,Sl')\0) (54) 



where U^x^Sl') = ^bUx, Si') - M\x,Sl'). Combin- 



ing (51) and (|54|) with (50) and performing the integra- 
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3/2 



lx(*)> = 



\fafn \ 

— — — e 2 (!+ 

1 + ihat/M j 



t/M) 11 — 



2M 



Vr-T) dt' -(R-Ro)- j dCl j dn 



OO pOO 



+ f £^ Xj Q} X tf^ a ^ d 3 x _ eo j as- [v(x,n,n')- ^t 3 Tr[v( x ,9.,n' 

dV 



)|0) (55) 



where, 



1 + ihat/M 




(56) 



The wave-packet evolves into a sum of two amplitudes: 
the first for the field remaining in the vacuum state and 
the centre of mass spreading as would be predicted in 
the absence of any coupling to the field; and the sec- 
ond for pairs of field excitations to be created, and the 
centre of mass to be shifted by a quantity equal to the 
double frequency integral over the oscillatory factor, ( 56 ) 



times a quadratic function of the expansion coefficients 
(see appendix |A|. The magnitude of this fluctuation in 
the wave packet evidently depends on the size, shape, 
and dispersion of the object. The relevant timescale of 
this fluctuation is one over the frequency where the imag- 
inary part of the susceptibilities is largest. Typically this 
timescale is many orders of magnitude faster than the 
spreading time of the wave-packet, T s . 



III. CONCLUSIONS 

It has been shown that the action ^ — previously de- 
rived in [38] — can be used as the basis for both a classical 
and quantum theory of radiation pressure, fully includ- 
ing the effects of dispersion and dissipation. The action 
is constructed so that it gives the Maxwell equations con- 
taining the correct constitutive relations for electromag- 
netism interacting with a moving medium. Here we in- 
vestigated varying this action with respect to the centre 
of mass coordinate to give equations of motion for the 
macroscopic body. 

We found that the expression for the classical force 
on a macroscopic body is given by (16 1. In the case 



where the body is stationary this reduces to the integral 
of the Maxwell stress tensor over the surface of the body, 
minus the time derivative of the Abraham momentum 
within the body. It was found that this is true even for 
dielectric bodies in contact, but that we must modify the 
action if we are to consider fluids or elastic media. 

Using the canonical procedure, a Hamiltonian was de- 
rived which was turned into a quantum mechanical op- 
erator and applied to the quantum motion of a small 
object. We found that the average force on the body is 



equal to the expectation value of the operator analogue 
of the classical result ( 45 ) . The detailed dynamics of an 



initially localised wave packet in vacuum were then anal- 
ysed. It was found that the peak of the wave packet 
is pushed by an amount equal to the usual expression 
for the Casimir force (in terms of electromagnetic Green 
functions), plus a time dependent factor corresponding 
to the creation of pairs of 'polariton' excitations within 
the body. The timescale of this fluctuation is determined 
by the frequency at which the object shows significant 
dissipation. 
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Appendix A: Normal mode expansion of field 
operators 

For reference here we give the expansion of the field 
operators that diagonalize ( 35 1 for a fixed position of the 
centre of mass R = Rq. This expansion is used in see- 
to diagonalize (35), giving (38). More details 



II A 



tion 

can be found on the derivation and verification of the va- 
lidity of this expansion within [321 1381 . A set of bosonic 
operators are defined with the following commutation re- 
lations, 



£i A) (*),c2'V) 



= t 3 d xy 5^(x-x')S(u}-u}') 



(Al) 

in order to define the number states associated with these 
operators, we write them in component form in terms of 
a set of unit vectors, e a 



cl X \x) = ^e (T CW(x) 



(A2) 



a set of states can then be defined in the usual way, 
with Ci x l(x)\0) = 0, cL X] J(x)\0) = \l u , a ^,x), etc. The 
canonical variables in ( 35 ) are expanded in terms of these 



bosonic operators. The expansion takes the following 
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form, 

A(x,t) 
where, 



dSl 



A(x, St)e~ iut + A' (x, St)e mt (A3) 



The curl operation given within O b is such that it takes 
derivatives with respect to the coordinate given in the 
argument. For example, 



O e (x, u) ■ G(x, x'; St) • O b (x', SI) = 



where the summation over A runs over the two indices, 
E and B (see section |l| . In the case of the reservoir of 
oscillators, the notation is only slightly more complicated 
than (A3-A4), taking the form, 



X^(x,t) 



where 



dSt 



X { ?(x,St)e 



nil 



X u (x,St)e 



lOI 



(A5) 



xf% X '>(x,x',oj,n).C"S'(x') 



(A6) 



In order that the change of variables ( A3 - A6 ) (along with 



the equivalent expressions for the canonical momentum) 
leave physical predictions unaltered, the canonical com- 
mutation relations, ( 32 33 1 must be preserved. One such 



expansion was found in |32j. which converts (J35J) into 
(38). In this case the expansion coefficients for the reser- 



voir operators are given by, 



(x, x',u, St) = USxx>S(u - St)5& (x - x') 

| d x (x,oj)-G(x,x';n)-d[,(x / ,n) 

(uj — St — irj) (lu + St + irj) 



(A7) 



where r\ is an infinitesimal quantity that serves to im- 
pose retarded boundary conditions on the operators. The 
canonical momenta operators of the reservoir are given 
by, 

fV£\x,x>,w,si) = -inf% x '\x,x',uj,n). (as) 



The two index quantity, G in ( A7 ) is the electromagnetic 



Green function, and is a solution of, 

V X ^(x,St)V X G(x,x';St)-St 2 e(x,St)G(x,x';St) 

= t 3 S^(x-x') (A9) 

where e = eo + Xee, and pT 1 — «n 1 — Xbb, with the sus- 
ceptibilities defined as in sectionhl The operators, O x (x) 
are given by, 



O e (x,uj) = iSta(x — R , w)1 3 
6 B (x,u)) = j3(x - R ,uj)Vx 



(A10) 



(A4) From (A10), the following identity can be derived which 



we give here for reference, 

J2o[(x,iiyd x (x,n) 



2st 



n a im[x M (x-jio,n)]i a 



xV -lm[x BB (x - R ,St)}V X 



(A12) 



Finally, the canonical variables of the electromagnetic 
field are expanded as, 



f$>(x, x', St) = G x (x, x'- St) • 0' x (x', St) (A13) 



and, 



x', St) = -iSte(x, St)G(x, x'; St) ■ 0' x (x\ SI) 

-l 3 a(x- Ro)5 xx ,S^(x-x') (A14) 



The subscript, '_L' in (A13) implies taking the transverse 



part with respect to the first index, and the first argu- 
ment, i.e. V • G±(x, x'; SI) = 0. 

For reference we also give the expansion coefficients for 
the electric and magnetic field operators defined below 
PI, 



f^(x,x',St) = iStG(x,x';St) • O x (x',St) 
f [x) (x,x',St) = V X G(x,x';St)-d[(x',St) (A15) 



using ( A15 ) we then find that the commutation relations 



between the Fourier components of the electric and mag- 
netic fields are, 



~t i tiSt 2 
5{x, Si),E (x', SI') = 6 (SI - St') Im[G(a;, x' , St)} 

J 7T 



and 



(A16) 



B(x,St),B (x',St') 
h 



S(St - Sl') h V x Im[G(x, x', St)} x & (A17) 



and finally, 



S(x,St),B\x',Sl')] = l —-5(St-St')Im{G(x,x',St)}xty 
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Appendix B: Regularization, and the vacuum force 
on a lone object 

Here we illustrate the canonical method for making 
the vacuum stress tensor finite, and demonstrate that 
the vacuum force on an isolated body in free space is 
zero. This vacuum stress tensor is given by 



He f°° ( 

lim / dQ Im< t(x, x' , fi) 



-~t 3 Ti[r{x,x',n)] \ (Bl) 



where t(x,x',Q) is defined by (52). As stated in the 
main text, (JbTJ) is not finite due to the behaviour of the 
integral over frequency. To remedy this, one usually sub- 
tracts Gq(x, x'\ fi) from the Green functions in (52), 



G(x, x'- fi) = G(x, x'; fi) - G (x, x'; fi) 



(B2) 



For piecewisc homogeneous media the quantity Go is the 
solution to (A9) for the case where e = const, and fi = 



const. The constant values of e and fj, are given by the 
local values where (Bl) is evaluated. This modification 
cannot affect the value of V • er , for lm[Go(x,x;fl)] 
does not depend on x (it is defined for a homogeneous 
medium) . 



We now take the imaginary part of ( Bl ) outside the in- 



tegral over frequency, and look to rotate the integration 
from the positive real axis to the imaginary one. This 
transformation is useful for the purpose of performing 
computations. As 101 — > oo in the upper half complex 



plane, it can be shown from ( A9 1 that G tends to a quan- 



tity proportional to a delta function in position, divided 
by frequency squared. This means that without apply- 
ing (B2 1 to (Bl ) both the integral over frequency and the 



limit of x — > x 1 diverge everywhere. However, subtract- 
ing Go from G removes this divergent term, leaving one 
with, at worst terms proportional to fi~ 4 in G. There- 



fore, once (B2) has been applied, the integrand of (Bl) 
goes to zero as |fi| -2 or faster in the upper half plane. 
This finding, along with the analyticity of the Green func- 
tions in the upper half plane mean we can thus rotate the 
integration from the positive real axis to the imaginary 
one, where the Green functions are real |35j . Taking the 
imaginary part of the resulting expression we obtain, 



er = hm 

7T X— <rx' 



t(x, x' , j£) 
1 



t 3 Ti[r(x,x',iO] 



(B3) 



where the bars indicate that ( B2 ) has been applied. This 



is the expression commonly used in calculations of the 
Casimir force [25]. 

We now show that when sufficiently isolated from other 
bodies, the force on an object computed from (B3) van- 



ishes. First it is shown that the divergence of (B3) van- 



ishes in free space. To see this we take the divergence in 



the following manner before the limit x — >• x' is taken, 
V-a = — lim / del(V + V) -t(x,x',i^ 

-i(V + V')'&[T(x,x / ) ie)]} (B4) 



To evaluate (B4) we note three things. Firstly, in free 



space G satisfies (A9) with /i = /iq and e = eo, and the 



right hand side set to zero. Therefore V-G(a;, x'; i£) = 0. 
Secondly, the gradient of the trace of rank two tensor 
equals, 



VTr[-r] = diTjj = e ijk E klm diT jm + dim 



= V ■ 



T X 



(B5) 



where repeated indices are summed over, and the 
star indicates that we have taken the dual of the 
two index quantity. Finally, due to reciprocity, i.e. 
G(x,x';Q) = G T (x , ,x-,n), we can interchange the two 
coordinates within the trace of r, VTr[r(a:, x' , ££)] = 
^Ty[t(x' ,x,i^)). Combining these three pieces of infor- 
mation, and dropping a term that vanishes in the x — > x' 
limit, 

V-d-o = 



lim 

TT x— >x' 



J cZ£ j f{x',x,iC) x V' 
t(x,x',i£) x ^ | 



(B6) 



Inserting the explicit form for r and applying ( |A9[ ) we 
can evaluate the dual quantities in (B6|, 



f(x,x',i£) x V' 



n 2 



(V-V')TT{G(x,x';iO] 



V -G(x',x;i£)- V ■ G(x,x';i£,)\ (B7) 



Inserting (B7) into (B6) and taking the limit x — ¥ x', we 



find the result that, in free space, 

V • &o = (B8) 
It was found in the main text that after inserting (TsTl) 



into ( 50 1 there is a term in the wave-function equal to 



the integral of V • <to over the dielectric body of interest. 
If the body is alone in free space (or at least very far from 
neighbouring bodies) then due to (B8) we can extend 
the integration domain over a very large region of space 
without altering the value of the integral. Therefore the 
value of the integral of V • CTq over the body is equal 
to J dS • <r over the surface of a sphere of large radius, 
centred on the body. To conclude we show that this 
vanishes. 

Writing G = G + Go and inserting this into (A9), we 
find that G obeys the differential equation 

V X fr 1 (x,iS,)V x G(x,x';ig) + € 2 e(x,i€)G(x,x'-,i£) 
= V X Xbb(zM0V X Ga(x,x';i£) 

-ex^{x,i^)G a {x,x';ii) (B9) 
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which has the solution 



G(x, x'; iO = J d 3 x"G(x, x"; if) • { x V" • Xbb(«)V" x -£ 2 x EE (i£)l 3 } ■ Go(*", a/; if) (BIO) 

I 



where Go(a:, £c'; fi) is the free space Green function. To 
leading order both G(x,x";fl) and Gq(x, x"; Q) go to 
••" _1 exp (— £|a; — a;"|) at large distances. 



zero as \x — x 



After the integral over £ in ( B3 1 has been performed one 
has a quantity going to zero faster than \x — x"\~ 2 , and 



J dS ■ &o therefore vanishes as the integration surface is 
taken to infinity. We thus conclude that, as anticipated, 
the integral of V • er over the volume of an isolated body 
is zero. 
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